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Abstract 

We discuss the behavior of a quantum glassy system coupled to a bath of 
quantum oscillators. We show that the system localizes in the absence of inter- 
actions when coupled to a subOhmic bath. When interactions are switched 
on localization disappears and the system undergoes a phase transition to- 
wards a glassy phase. We show that the position of the critical line separating 
the disordered and the ordered phases strongly depends on the coupling to 
the bath. For a given type of bath, the ordered glassy phase is favored by a 
stronger coupling. Ohmic, subOhmic and superOhmic baths lead to different 
transition lines. We draw our conclusions from the analysis of the partition 
function using the replicated imaginary-time formalism and from the study 
of the real-time dynamics of the coupled system using the Schwinger-Keldysh 
closed time-path formalism. 



1 



1 Introduction 



The effects of a dissipative environment on the dynamics of quantum systems have 
been intensively investigated during the last two decades ^. The most widely 
studied problem is that of a single macroscopic variable coupled to a set of micro- 
scopic degrees of freedom that act as a bath. The environment is usually described 
in terms of its collective excitations (lattice vibrations, spin or charge fluctuations, 
etc) that may be thought of as an ensemble of independent quantum harmonic os- 
cillators 0, |, ^ ^, 1^. Their coupling to the system is given in terms of a spectral 
density oc a cu* for <^ Uc, where a is a dimensionless coupling constant and 
iUc a high frequency cutoff. The exponent s characterizes different types of environ- 
ment. The Ohmic case (s = 1) is quite generally encountered but superOhmic 
(s > 1) and subOhmic (s < 1) baths also occur e.g. in the case of the Kondo effect 
in unconventional hosts 0, ^. 

The question of how dissipation destroys quantum coherence [|l], ^, |^ in two-level 
systems (tls) has been extensively investigated in the literature. The low-energy 
physics of many tunneling systems is well described by the spin-boson model §]. 
In this model, the two equivalent degenerate states of the tls are represented by the 
two eigenstates = ±1 of an Ising pseudo-spin. A transverse field coupled to 
(say) represents the tunneling matrix element. Much is known about the properties 
of this model and its relationship to several other models including the 1-D Ising 
model with inverse squared interactions |jlO|], the anisotropic Kondo model 0, |l^ 



or the resonant model |T^. Three different regimes are possible depending on the 
value of a: in the Ohmic case, at zero temperature, there is a phase transition at 
a = 1 For a < 1 there is tunneling and two distinct regimes develop. If 

a < 1/2 the system relaxes with damped coherent oscillations; in the intermediate 
region 1/2 < a < 1 the system relaxes incoherently. For a > 1 quantum tunneling 
is suppressed and (cr^) 7^ signalling that the system remains localized in the state 
in which it was prepared. 

These results also hold for sub-Ohmic baths while weakly damped oscillations 
persist for super-Ohmic baths 0] . At finite temperatures (but low enough such that 
thermal activation can be neglected), there is no localization but the probability 
of finding the system in the state it was prepared decreases slowly with time for 



These conclusions, derived for a single tls interacting with a bath, can be applied 
to a macroscopic system in the diluted regime, i.e. when the interactions between 



the TLS are unimportant compared with those between a tls and the bath |16]. 
There are, however, physical systems that can be viewed as a dense set of tls in 
which their mutual interactions can no longer be neglected. The question then arises 
as to which are the effects of the interplay between the interactions between the tls 
and their coupling to the noise on the physics of the interacting system. 

In this paper we discuss this issue in the context of a glassy macroscopic system 
with random, long-ranged interactions. This situation is realized experimentally in 
systems such as uniaxial spin glasses in a transverse magnetic field and disor- 



dered Kondo alloys [T3, ITUl . Metallic glasses with tunneling defects are also systems 
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in which the effects that are of interest here could be observed experimentally. 

In thermodynamic equilibrium, in the absence of the bath, the interactions be- 
tween the TLS lead to the appearance of an ordered state at low enough temperature. 
If the interactions are of random sign, as in the models we consider here, the lat- 
ter will be a spin glass (sg) state. In this phase the symmetry between the states 

= ±1 at any particular site is broken but there is no global magnetization, 
J2i{<^i) = 0. Since the coupling to the bath also tends to locally break the symme- 
try between the degenerate states of the tls, both interactions compete with the 
tunneling term in the Hamiltonian. We thus expect the presence of noise to increase 
the stabihty of the SG state against quantum fluctuations. The consequences of 
this fact are particularly interesting when the coupling to the bath leads by itself 
to localization at some a = a*^^". Consider a system of size with a > a*^^'^ at 
T = and suppose that we turn off the interactions between the tls. The ground 
state of the system is then 2^-fold degenerate as each tls can be in one of the 
states (erf) = ± (To (say) independently. If we now turn on an infinitesimal random 
interaction between the tls, this macroscopic degeneracy will be immediately lifted 
as the system will select among its 2^ degenerate configurations the one (or one 
among the ones) that minimizes the interaction energy. If we denote by J the typi- 
cal scale of the interactions and by a'^''" the localization threshold, we thus expect 
a quantum critical point at J = 0, a = between a quantum paramagnet and 
the ordered state such that, for a > 0;*^^", the SG phase survives down to J = 0. 

A system of non-interacting localized tls and a SG state in equilibrium are 
in some way similar: in both cases ^i{(y't) = and the presence of order is re- 
flected by a non-vanishing value of the long-time limit of the correlation function, 

= ^^^t^oo J2i{(^iit)'^i{^)) (since we assume equilibration the correlation is 
stationary and the reference time can be taken to be zero). However, this resem- 
blance is only superficial. In the renormalization-group language, J is a relevant 



variable |]20|. Therefore the details of the dynamics of the two systems are expected 
to be quite different, in particular the way in which the correlation function reaches 
its asymptotic limit, ^ea, that determines the low-energy part of the excitation spec- 
trum of the system. 

Further differences between the localized state and the SG state are seen from 
the study of the out of equilibrium relaxation of such states. Indeed, an impor- 
tant feature of glassy systems is that their low-temperature dynamics occurs out of 
equilibrium. If the system is macroscopic, its size is very large (diverges in the 
thermodynamic limit). In a realistic macroscopic situation, the asymptotic long- 
time limit follows this large size limit. Many experiments, simulations and analyti- 
cal studies show that the time needed to reach equilibrium after entering the glassy 
phase diverges so quickly that the relevant relaxation occurs out of equilibrium. The 
dynamics at low temperatures is then non stationary, i.e. the dynamic correlation 
functions loose time translation invariance. If denotes the time elapsed since a 
quench from the high temperature phase into the SG phase, C(t + 1^,1^) depends 
on both t and ty^. The order in which the limits ^ 00 and t 00 are taken is in 
this case very important. For sufficiently long t and but in the regime t <^tyj, the 
dynamics is stationary and the correlation function reaches a plateau ^ea- Much of 
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what was said above for the equihbrium state also holds for this stationary regime. 
However, for times t > t^, the system enters an aging regime where the correlation 
function depends on the waiting-time t^, explicitly. In this regime, the dynamic 
correlation function vanishes at long times, lim^^oo C(t + t^, ty^) = 0, at a rate that 
depends on tw (This is to be confronted to the dynamics in the localized state, 
where C(t + tu,,tu,) reaches, for any waiting-time and long enough t a plateau 
that it never leaves.) In this regime, even for a > a'^'''*, small interactions will result 
in the destruction of localization of the tls at long enough times. 

The problem of a single tls being a difficult one, that of an infinite set of interact- 
ing TLS seems hardly solvable at this stage. Therefore, as a first step, we shall focus 
on the low-temperature dynamics of a very simple model that mimics some of the 
features of more realistic ones. This is a quantum generalization PTI E^, 123, E3, El] 



of the random p-spin spherical model I^Bj coupled to a bath of quantum harmonic 
oscillators. The principal merit of this model is that it is simple enough that it can 
be studied in detail. Yet, many of its properties are generic and expected to hold at 



least qualitatively for more realistic models [27|. The usual methods of equilibrium 



quantum statistical mechanics are inappropriate to describe the nonstationary situ- 
ation. We solve the model using two methods especially designed to treat systems 
out of equilibrium. 

One is based on the Schwinger-Keldysh real-time approach to non- stationary sys- 
tems. It was first applied in this context to the quantum p-spin model in Ref. 
and used subsequently in other cases including the SU(A/') fully connected Heisen- 
berg model in the limit of large J\f and the soft spin version of the Sherrington 
- Kirkpatrick model [^. It allows one to obtain the full time dependence of the 
correlation C(t + t^, t^). 

The second method is based on the Ansatz of marginal stability (ams) within 
the replica analysis of the partition function. Originally developed for classical sys- 
tems [0, this method was recently used to discuss the low- temperature properties 
of quantum glassy systems [0, |3^. Its main advantage is that it uses a formal- 
ism that is closely related to the imaginary-time approach to equilibrium quantum 
statistical mechanics. In Ref. |2i| the ams was extensively applied to the quantum 
spherical p-spin model in the absence of the bath. It was shown that the position of 
the dynamic transition line predicted by this method coincides precisely with that 
obtained using the real-time approach. It was also shown that the time dependent 
correlation function computed using the ams in the absence of the bath is identical 
to the stationary part of the non-equilibrium correlation function (C > ^ea) when 
one takes the long-time limit first and the limit in which the coupling to the bath 
goes to zero next. The marginality condition imposed by the Ansatz is intimately 
related to the fact that the correlation will further decay from ^ea towards zero. 
(The details of this second decay as, for instance, the two-time scaling are not ac- 
cessible with this method.) A localized solution with C(t-t-t^,t^) approaching, and 
never leaving, the plateau at gEA corresponds, in replica terms, to a stable replica 
symmetry solution. In this paper we extend the ams to study the dynamics of the 
model in the case in which the system is coupled to the environment. 

This paper is organized as follows. In Section 1^ we motivate and introduce our 
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model and discuss its relationship to the more usual spin-boson model. In Section 3 
we outline the imaginary-time formalism used to solve the problem at equilibrium 
and within the AMS. We compute the partition function of the coupled system 
and determine its phase diagram in both situations. We also discuss the long-time 
dynamics of the coupled system using a very accurate long-time approximation 
p5[] that allows us to solve the model analytically. This approximation is then used 
to discuss the influence of a coupling to different types of environment on the T = 
quantum phase transition. The real-time dynamic of the system is discussed in 
Section H and the results are compared with those of Section 3. Section ^ contains 
a brief summary of our main results and our concluding remarks. 



2 The model 



In order to motivate our model, we start by considering a collection of identical 
interacting TLS coupled to a bath of independent harmonic oscillators P, [T^. We 
assume for the moment and until otherwise stated that the combined system is in 
thermodynamic equilibrium. The Hamiltonian of the coupled system may be written 
as 

H = Hs + Hb + Hsb, (2.1) 

where Hs, Hb, and Hsb denote the Hamiltonians of the system, the bath and their 
coupling, respectively. These are given by: 



N 



Hs 

Hb 
Hsb 



- H 4 xe^! ■ 



(2.2) 

(2.3) 
(2.4) 



i.e 



Here, the Pauli-matrices erf represent the tls's pseudo-spins, A/h is their tunneling 
frequency and V their mutual interaction potential that we leave unspecified for the 
moment. Xi and are the coordinate and momentum of the £-th oscillator and me 
and ue its mass and frequency, respectively. We denote by c} the coupling constant 
between the i-th tls and the i-th oscillator. 

the oscillator degrees of freedom may be inte- 



Using standard methods |]J 
grated out to express the partition function of the system solely in terms of the tls 
variables as: 

with 



r exp ( -- 



(2.5) 



S = l dr}^-Aj:^nr) + V[a\r)]^ + -j: I dr dr' K^.ir-r') a,\T)afiT') , 

(2.6) 
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where T the imaginary-time ordering operator and we have introduced the notation 
The kernel Kijir) in Eq. is |14| 

Kijir) = Kijiuk) exp(-iu;fcr) , (2.7) 

where ojk = 2-nk/{h(3) are the Matsubara frequencies, the coefficients Kij{ijjk) are 
given by 

A.(,,).r^M£)_^. (2,8) 

and we have introduced the spectral density of the environment hjiuj) through 

hAuo) = TT V ^^^^l±^J^5{u - UJ,) . (2.9) 

We make the simplifying assumption that the dynamic interaction between different 
TLS generated by integration over the degrees of freedom of the bath can be neglected 
compared to the static interaction potential included in ^[c?^]. Therefore we write 

Ii,{uj) = 6^,I{uj) , (2.10) 

and we choose the standard parametrization 



liuj) = 2ah{ — \ CO e"'^/"^ , (2.11) 

where a is a dimensionless coupling constant, Uc is a high frequency cutoff and 
is a microscopic phonon frequency necessary in the non-Ohmic cases in order to 
keep a dimensionless. For simplicity, we shall restrict the exponent s to lie in the 
interval < s < 2. With this choice the integral on the right-hand side of Eq. (|2.8| ) 
converges without the need of introducing an infrared cutoff and the upper cutoff 
may be eliminated by taking the limit a;c — >^ oo. This leads to the expression: 

KM = ."^^ M' . (2.12) 

^ph sm(7rs/2) 

We shall consider diagonal p-spin interactions of the form 

N 

V[a']= E Jn...,<-K^ (2.13) 
n<...<jp 

with random couplings Ji^...ij,- These are taken from a Gaussian distribution with 
zero mean and variance 

^n...v)' = ^~V/(2iVP-^), (2.14) 
where the overline represents an average over disorder. 



s-l 
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In the case p = 2 and for an Ohmic bath (s = 1), Eq. ( |2.6| ) is equivalent to 
the action of a disordered Kondo alloy model ||33|. In this context ||12[, A = J^, 
the transverse Kondo coupling and (1 — a) <^ 1 is proportional to Jy^ the parallel 
Kondo coupling. In the opposite limit, a <C 1, Eq. (p.6|) is a representation of the 



partition function of the SK spin-glass model in a transverse magnetic field 
weakly coupled to a phonon (or spin [0) bath. 

The main difficulty in solving the quantum statistical mechanical problem defined 
by Eq. ( |2.6| ) stems from the discrete nature of the spins. It was shown in Ref. 
(hereafter referred to as CGS) that, in the absence of the bath, a solvable (and yet 
non-trivial) model can be obtained by generalizing the eigenvalues Sj = ±1 to 
continuous variables — oo < Sj < oo, and replacing the hard constraint = 1 by 
the soft spherical constraint J^iisf) = N. The derivation of the effective continuous 
model in the presence of the bath is analogous to that given in CGS for the isolated 
system and we refer the reader to this reference for the details. After performing 
a Trotter-like decomposition of the time-ordered exponential in Eq. ( |2.5| ) the trace 
becomes a functional integral over classical fields Sj(r) and the tunneling term in 
the action acquires the low-energy form 

l\/r ^KR / a„z \ 2 

(2.15) 




where we introduced the mass 

and To is a cutoff representing a microscopic spin-flip time that we identify with 0J~^. 
Since we work in the regime in which hujc is the highest energy scale in the problem, 

< M < cx). 

The continuous version of Eq. ( p.6|) is thus given by 



1 r^P 



^ -0 



M + fdr'Kir - r') s,{r)s,{r') + '-{s]{r) - 1) 



- H dr Ji^,„i^Si,{T) . . .Si^ir) , (2.17) 

h<...<ip "^^ 

where 2; is a Lagrange multiplier that enforces the spherical constraint 

{s{t) . s{t')) |_, =^J:{\^M\') = N , (2.18) 
lip 



where the angular brackets represent the average with respect to the action (|2.17|) 



Equations (p.l7|) and (|2.18|) deflne the quantum p-spin spherical model that we 



discuss in the rest of this paper. The mass parameter M is a measure of the strength 
of quantum tunneling. If Atq/H ^ 1, M is large. In this case, the gradient term 
favors conflgurations in which Sj(r) is almost r- independent. The partition function 
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is then largely dominated by the contribution from the static fluctuations {i.e. those 
with = 0). Since K{ujk = 0) = 0, these variables are unaffected by the coupling 
to the bath which drops out of the partition function in this limit. With increasing 
A, M decreases and the amplitude of the quantum fluctuations becomes large. The 
r-dependence of Si(r) then becomes essential. 

There are three points worth discussing before presenting the solution of the 
model. The first one is their dependence on the value of p. For p = 2 the action is 
quadratic and the problem is readily diagonalizable by Fourier transformation. This 
simple case was extensively discussed in the literature both without [35| and with [3^ 



a bath. In the former case, the competition between the mass and interaction 
terms in Eq. (|2.17|) leads to the existence of a critical mass ~ ft /J above 



which the ground state of the system overcomes quantum fluctuations and acquires 
glassy order. For p = 2, however, the ordered phase is of a trivial type, with a 
structureless order parameter (see below). Its physical properties are non-generic 
and qualitatively different from those of discrete spin systems. The presence of a 
coupling to the bath does not change this situation. For all p > 2 the ordered ground 



state is non-trivial |2^, ^] and the model shares a number of qualitative 

features with more realistic ones. Therefore, from here on, we shall discuss this case, 
choosing the particular value p = 3 in our numerical calculations. 

The second point is about the case J = 0. In this case Eq. ( p.l7| ) reduces to 
a simplified model for a tls whose physics differs in some ways from that of real 
two-level systems. For J = we have 

{\sM\')j.,= , ^ , (2.19) 



where K{uJk) is defined in Eq. ( |2:T2| ) and Eq. ( p38|) at T = reads 



1 /""^ Ti 

1 = - / duj = = fsiz) . (2.20) 

Trio Muj^ + K{uj) + z ^ ^ ^ ^ 

We consider the Ohmic case first. For s = 1, Eq. (|2.19|) is the propagator of a simple 
damped harmonic oscillator with frequency c^o = \J z/M self-consistently determined 
by Eq. (p.20| ). From the position of the poles of Eq. ( |2.19|) we see that there is a 



transition between underdamped and overdamped regimes at z = cP'Ti^ j (4M). Using 
this value of z in Eq. (|2.20| ) (with s = 1) and solving for a we find that this occurs 
at a = 2 /it, independent of M. Away from this value we easily find the following 
limiting behaviors: 

(2.21) 

2 2 

exp (— vra) , a 1 

For a -C 1 the system exhibits weakly damped oscillations with frequency ujq ~ 
h/M. In the opposite limit, a ^ 1, the correlation function decays exponentially 
with a time-constant that increases exponentially with the strength of the coupling, 
r ~ exp (vra). 
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Comparing with the resuhs for the spin-boson model summarized in the Intro- 
duction, we see that the transition between coherent and incoherent motion at a 
universal value of a < 1 is preserved in the spherical model but the localization 
transition at a = 1 is replaced by a crossover to a high coupling regime character- 
ized by an exponentially small energy scale cx exp (— vro;). In this regime, tunneling 
is not suppressed but its rate is strongly reduced. 

In the superOhmic case there is no localization transition either. One can easily 
show that for s > 1 the expression on the second line of Eq. ( |2.21| ) is replaced by 
z ~ ^cUp/jO^/*-^"'*'' for a ^ 1. The decay rate of the correlation function still decreases 
continuously as the strength of the coupling to the bath increases but only as a power 
law. 

In the case of a subOhmic environment the situation is different. For s < 1 
the integral on the right-hand side of Eq. ( |2.20| ) is finite at z = where it takes 
its maximum value. For /s(0) < 1 Eq. ( p. 20 ) cannot be satisfied for any positive 
value of z. This phenomenon, completely analogous to Bose-Einstein condensation, 
signals a localization transition. Solving the equation /s(0) = 1 we find the critical 
coupling given by 

~ {^^—\ ■ (2-22) 

Conversely, for any value of a the system localizes for a sufficiently high value of M. 
These results are analogous to those obtained for discrete tls 

One must keep these differences between the original model and its spherical 
version in mind when interpreting our results, especially those discussed in Sec- 
tion p.4.2| . In Section 4 we shall illustrate the interplay between the localization 



observed for J = and and the glassy dynamics that appears when J > 0. 

The third point we want to stress is that the coupled system can be thought 
as describing the motion of a quantum Brownian particle of mass M, constrained 
to move on a A^-hypersphere of radius \/N , in the presence of a random potential 
V{s). The Brownian nature of the motion arises because of its interaction with 
the quantum thermal bath. The infinite dimensional spherical limit yields, however, 
unphysical results if one wants to compare it to the well-known problem of the diffu- 
sion of a free quantum particle coupled to a phonon bath in a finite D dimensional 
space. While here we find a localization transition when the bath is subOhmic, 
such a transition does not exist in the absence of disorder in the finite dimensional 
problem. 



3 Replica solution 

It is by now well established that several properties of disordered systems can be 
derived with the help of the replica trick. This approach enables one to derive an 
effective action for an imaginary time dependent matrix order parameter. It has 
been noticed that different Ansatze that parametrize this order parameter describe 
different physical situations as thermal equilibrium (equilibrium condition) and the 
asymptotic dynamic regime {Ansatz of marginal stability - AMs). The bulk of this 
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Section is devoted to the analysis of the consequences of the AMS. The definition of 
this Ansatz excludes localization as a possibility. As discussed in the Introduction 
we expect that as soon as interactions are switched on, full localization is replaced 
by a glassy solution with non-trivial dynamics. This argument justifies the use of the 
AMS from the start. We briefly comment at the end on the equilibrium properties 
of the model. 



3.1 Formalism 

The presence of disorder makes it necessary to compute the averages of all physical 
quantities and, in particular, of the free-energy. To this effect we use the replica 
trick, i.e. we write 

1 



/?/ = -— InZ = -— lim-lnZ" . (3.1.1) 

The derivation of the expression for the free-energy associated to the action in 
Eq. ( |2.17| ) closely follows that performed for the isolated system in CGS where the 
interested reader will find all the necessary details. An imaginary-time dependent 
order parameter Qab{T, t') is defined as 



Qab{r,T') = -{Ur)-s,{T')) , (3.1.2) 

where a, b are replica indices. The spherical constraint imposes the restriction 
Qaa(O) = 1. We are interested in a stationary situation in which (5ab(T, t') depends 
only on time differences and is a periodic function of its argument with period (3h. 
We thus introduce the Fourier transforms Qabi^^k) = /o^^ drQabiT) exp(2ti;fc) in terms 
of which the averaged free energy is found as 



where 



/5/=limGo, (3.1.3) 



2^0 = --ETrln[(/3^)-^Ql -5:fl-4-E^a6(^fc)Qa6(^fc) 



k k \ ab 



■^EX dr[^-^eM-^^ur)Q^^M) - , (3.1.4) 



and the operator Oab{,^k) is defined by 

OabM = ~l5ab {MujI + Z + k{uJk)) ■ (3.1.5) 

The equations of motion are found from the saddle-point of the free-energy with 
respect to variations of Qabi^^k)- They read: 

-{Mul + z + K{uju))6ab = (Q-^)^^ M + dTexp{zujkT)Ql,\r) . (3.1.6) 
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Equation (|3.1.6| ) together with the spherical constraint QaaiT = 0) = 1 determine 
the different phases in the modeL 

In the following, we discuss the solutions to Eq. ( p.l.6|) . Except when otherwise 
stated, we shall work with dimensionless variables. These are defined by measuring 
energies in units of J and time in units of h/J. The strength of quantum tunneling 
and of the coupling to the bath are then measured by the parameters 



a 



{JM) 



a. 



( 



ph 



sm 



(vrs/2)V J 



respectively. 

In the paramagnetic phase (pm) Qabi^k) is a diagonal matrix, 

Qab{^k) = qd{^k)5ab ■ 



(3.1.7) 



(3.1.8) 



Equation ( p.l.6| ) then reduces to 



^k 



1 



2 



(3.1.9) 



In the SG phase, we search for 1-step replica symmetry breaking (rsb) solutions of 
the form 

Qabi^^k) = {qdi^^k) - Qea) + qEA^ab, (3.1.10) 

where eab = 1 if a and b belong to the same diagonal block of size m x m and 
zero otherwise, and we introduced the Edwards- Anderson order parameter ^ea- It 
was shown in CGS that this Ansatz is an exact solution of the isolated model. The 
proof still holds in the presence of the bath provided that Xim^^^Q K {uj) = 0, which 
is verified here (cf. Eq. (|2.12|) ). 

To completely determine the order-parameter matrix, ggA and m must be com- 
puted. As discussed in detail in CGS, this may be done in two different ways, each 
leading to a physically different state. Within the Ansatz of marginal stability, q^A 
is determined by extremization of the free energy, m is chosen such that the sta- 
bility of the ordered state is marginal, i.e., that its excitation spectrum contains a 
zero-energy mode. 

Decomposing the diagonal order-parameter qd{uj) in a singular and a regular 
part. 



qd[^^k) 



(3.1.11) 



an equation for gREG('^fc) can be derived by a straightforward generalization of the 
results of CGS to the case in which noise is present. It reads: 



(tUfc) - Sreg(0 



5'reg('^A:J 



with 



z = -prnqP^ — 



(3.1.12) 



(3.1.13) 
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(3m 



Qea 



and 



Xp 

m + Xp 



-p/2 



dr (g^-i(r) 



"iEA 



COS [UJkTj 



The parameter Xp takes the value (see CGS for the details of the derivation) 



Xn 



V 



2 . 



(3.1.14) 
(3.1.15) 

(3.1.16) 

(3.1.17) 



For given p and a, Eqs. (|3.1.12| )-( |3.1.16| ) have solutions with q^^^ ^ only for low 
enough values of F and T. Otherwise, thermal or quantum fluctuations destroy the 
ordered state. There is thus a boundary Fc(T) in the T — F above which the system 
is in the PM state. We determine its shape below. 



3.2 The dynamic phase diagram 

We determined the phase diagram for the coupled system for p = 3 using the 
numerical methods described in CGS. A critical line with a second order section 
(close to the classical critical point (T^^, F = 0)) and a first order section (close 
to the quantum critical point (T = 0,Frf)) is also obtained in the presence of an 
environment. The second order critical line is determined by the condition m = 1, 
the first order critical line is defined as the locus of the points where a marginally 
stable solution first appears with decreasing F for T fixed (see Fig. |^). For each 
F and a this defines a dynamic transition temperature Trf(F,a). It was shown in 
CGS that r^(F, a) precisely coincides with the temperature below which the real- 
time dynamics of the system looses time-translation invariance and the fiuctuation- 
dissipation theorem (fdt) is violated ||22| . 



The qualitative features of the phase diagram, similar to those found for the 
isolated system, are as follows. For p > 2, the transition is discontinuous in the 
sense that the order parameter ^ea jumps across the phase boundary. The transition 
line contains a tricritical point (T*, F*) that divides it in two sections. For T > T*, 
physical properties are continuous across the transition. The latter is therefore 
second order in the thermodynamic sense. For T < T*, instead, physical quantities 
are discontinuous across the transition which is thus first order. The origin of this 
behavior is the fact that the values taken by the parameter m on the transition line 
are different above and below T*. For T > T*, m = 1 along the transition line. 
This is its value in the paramagnetic phase meaning that m is continuous across 
the transition and so are the observables. For T < T*, m ^ 1 along the transition 
line but it is a decreasing function of T that vanishes linearly as T — 0. Crossing 
the phase boundary at T < T*, m is discontinuous and so are physical properties. 
Notice that the line Trf(F,a) lies always above Ts(F,a), the static critical line that 
we shall discuss below. 

We show on the right panel of Fig. |l| the dynamic phase diagrams obtained for 
p = 3 and three values of the coupling to an Ohmic bath, a = 0, 0.25, 0.5. The full 
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Figure 1: Static (left) and dynamic (right) phase diagrams for the p = 3 spin model 
coupled to an Ohmic bath (s = 1). The couplings to the bath are a = 0, 0.25, and 
0.5 from bottom to top. The solid line and line-points represent second and first 
order transitions, respectively. 

line and the line-points represent second and first order transition, respectively. We 
make the following observations: 

1. In the limit F ^ the transition temperature is independent of the strength 
of the coupling to the bath. 

2. The size of the region in phase space where the system is in the ordered state 
increases with a. Coupling to the dissipative environment thus stabilizes this 
state. 

3. The dynamic tricritical temperature decreases rapidly with increasing a. 

Our first observation is a consequence of the fact that in the limit F — the 
partition function is essentially determined by the zero-frequency components of 
the pseudo-spin which are decoupled from the bath (see Section 2). This result is 
however non-trivial from a dynamical point of view, since it implies that the dynamic 
transition of a classical system coupled to a colored classical bath is not modified 
by the latter. 

The second feature follows from simple physical considerations. The interaction 
term in the action favors spin-glass order. Coupling to the bath favors localization 
and its effect is to reduce the effective tunneling frequency. Therefore, in the presence 
of the bath, the value of the bare tunneling frequency needed to destroy the ordered 
state must increase with a. Even if the localized state and the glassy state may seem 
superficially similar, their are indeed very different. In the former, the correlation 
function C(t + t^^ty^) approaches a plateau as a function of t and never decays 
towards zero while in the latter the relaxation first approaches a plateau but it 
eventually leaves it to reach zero for t 3> We shall see this difference explicitly 
in the analysis of the real-time dynamics of Section 4. 
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Figure 2: The diagonal part qd{T) for the static (left) and dynamic (right) solutions 

The fact that the coupling to the environment favors the ordered state also 
reflects itself in the value taken by the order parameters ^'^(t) and ^ea- We display 
in Fig 1^ the r dependence of the diagonal part of the order parameter, ^^(t) for the 
static and dynamic solutions at a fixed temperature and F for different values of 
a. It can be seen that, as a increases, ^^(t) reaches a higher plateau level at long 
imaginary times. The analysis of q^^ is postponed to Section 3.4 (see Fig. 5). 

Figure | displays the m-dependence of F at a fixed temperature (T < T*), for 
different values of the coupling to the noise. The function F(m) is double-valued 
and the physical branch is that on which dm/dV > 0. This is a consequence of 
Eq. ( |3.1.13| ) that shows that m is a decreasing function of ^ea which itself is a 



decreasing function of F. It can be seen that for fixed F and T, m decreases with 
increasing a. Thus, the coupling to the bath results in a higher effective temperature 
in the glassy phase (see Section 4 for a definition of Teff and a discussion on this 
issue). 

We have also studied the phase diagram in the non-Ohmic cases. Figure ^ shows 
a comparison of the effects of an Ohmic bath and two non-Ohmic ones, subOhmic 
{s = 1/2) and superOhmic (s = 3/2) for the same value of a. It may be seen that 
for the chosen values of the parameters the region of stability of the ordered phase 
is enhanced (reduced) for a subOhmic (superOhmic) bath with respect to an Ohmic 
one. This feature is not generic as there are other values of for which the relative 
sizes of the effects of Ohmic and nonOhmic baths are different. Indeed, in preparing 
these figures we used ujph = 10 in the nonOhmic cases and this parameter modifies 



the coupling to the bath due to the factor cj^ ^ in I (to) 



3.3 Equilibrium phase diagram. 

In thermodynamic equilibrium both ggA and m are determined by imposing that 
the free energy is an extremum with respect to their variation. This leads to the 
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Figure 3: The break point m as a function of F for three values of the couphng to 
an Ohmic environment a. (3 = 20 > (3*. 



conventional thermodynamic equilibrium state. The value of Xp is obtained from, 

2 

- log(l + Xp)+ .^^^ =0. (3.3.1) 

The transition line is defined as the locus of the points where the free energies of the 
PM and SG phases coincide. For each F and a this defines a freezing temperature 
Ts[r,a) at which the system enters de SG state. The qualitative features of the 
equilibrium phase diagram shown on the left panel of Fig. ^ are similar to those 
found for the dynamic case. Notice that the line Trf(F, a) lies always above Ts(F, a) 
and that, in contrast to what we found for the dynamic tricritical temperature, the 
equilibrium tricritical temperature T* depends only weakly on the strength of the 
coupling to the bath. 



3.4 Low-energy properties of the marginal SG state 

Insights on the low-energy properties of the model may be gained by studying it 
in the framework of a simple and accurate approximation applied to the isolated 
model in CGS. It consists in deriving the exact low- frequency form of qREci^k) and 
using it over the whole frequency range assuming that physical properties at low 
temperatures are mainly determined by the low-energy excitations of the system. 
We consider in the following the T = case. 

3.4.1 The low-frequency form of q-REei^k) 

We start by assuming (and verifying later) that qREcX'^) (cf- Eq. ( p.l.ll|) ) decays in 
imaginary time as a power-law: 

qKEcXr) « |r|-^ . (3.4.1) 
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Figure 4: The dynamic critical line for subOhmic {s = 1/2, upper curve), Ohmic 
(s = 1, middle curve) and superOhmic (s = 1.5, lower curve) baths. huph/J = 10 
in the non-Ohmic cases. The dimensionless coupling to the bath is a = 0.5 in all 
cases. 

Then, we may write (cf. Eq. ( ^.1.16D ): 



^HEG (tUfc) - Sreg(O) 



= ^ Tdr (cosa;,r-l)[(p-l)g^;VEG(r) + ...|3.4.2) 
2 Jo L -f 

oc \uJkf~^ (IH ) , 



where the dots represent terms that vanish in the limit Uk 0. Therefore, in the 
long-time limit, 

t^^aM - SaEG(O) ^ [gREG(^) - gREG(O)] (3.4.3) 



Substituting Eq. (|3.4.3|) in Eq. (|3.1.12| ) and solving for Qregi^) we find: 

2 



ul/T + a^liVkl-' + 2/tp + ^Jul/T + as\uJk\yujl/T + a^kfel* + inp 

(3.4.4) 

where we introduced the parameter 



Kp — 



lp{p — l)q 



p-2 
EA 



(3.4.5) 



Equation ( p.4.4|) only holds in the low-frequency limit where it reduces to 



1 - 4= + 



1/2 



-1 



/«s I is/2 



(3.4.6) 
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This leads to the long-r behavior 



qKEcXr) ~ . (3.4.7) 



The assumption ( |3.4.1| ) is thus self-consistent with the exponent ( = 1 + s/2. In the 



absence of the bath Eq. (|3.4.4| ) leads to 



QREoir) ~ -j= , (3.4.8) 



the result found previously in CGS for the isolated system. A crossover between 
these two regimes occurs at Tcr = (Ftts)^^^'*"^''. The analytic continuation of Eq. ( p.4.6| ) 
yields the low-frequency limit of the imaginary part of the susceptibility: 



1/2 
1 1/2 



xLo(^)~sign(c.)(^) |^r/^ (3.4.9) 



The result for the Ohmic case, Xwai^) oc Ic^l ' was previously found for p = 2 
continuous and discrete Kondo-alloy models at the quantum critical point |3^. In 
the marginally stable state, this behavior persists throughout the low-temperature 
phase. 

Equation ( p.4.4| ), exact in the limit Uk —>■ 0, may be used as an approximation 
for finite frequencies. It will be seen in the following that this approximation allows 
one to gain useful insight on the effects of the environment on the physics of the 
interacting system. 

3.4.2 The quantum phase transition 



The normalization condition and Eq. ( 3.1.11|) lead to the following equation for the 
order parameter at T = 0: 

1 - ?EA = ^ 51 ^REG(^fe) gREG(^) , (3.4.10) 



This is still an implicit equation for the order parameter as gREG(i-^) depends on 
through (cf. Eq. (|A5D). 



We now approximate Eq. ( |3.4.10| ) by assuming that the integral on the right-hand 
side is dominated by the low frequencies. Then, we use for gaEG(i^) the expression 
given in Eq. ( p.4.4| ) which is asymptotically exact in this limit and write: 



9 fli,i 
l-q.. = - , , (3.4.11) 



It is convenient to make the change of variables oj = (TasY^^'^ Eq. ( |3.4.11| ) 

which leads to: 

^« (1 - = / 2^ s^n ^ nr^s / 2^ sr T = 3s{e) , (3.4.12) 
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where 



Equation ( |3.4.12|) will be used to study the T = quantum phase transition. We 
shall mostly be interested in the vicinity of the quantum transition where the system 
is close to the quantum paramagnetic state. We discuss separately different types 
of environment. 



a 



Mu. 



ph 



sin^ ns/2 



(3.4.14) 



The Ohmic case 



Setting s = 1 in Eq. ( |3.4.12| ) the equation of state may be written as: 



—(1 - = / nr^-^ T = ^1^')' ^^-^-^^^ 

z JQ + X + 2e + \/x^ + XV + X + 4e 

with e = Kpl (Fq;^). 

We show in Fig. |^ the a-dependence of ^ea in the marginally stable case for 
p = 3 at fixed F = 4. We represent with line-points the results obtained numerically 
from the full equations at a finite but low temperature, T = 0.1. The dashed 
line instead represents the approximate solution derived from Eq. (|3.4.11|) . The 



agreement between the two calculations is very good even if the approximation 
strictly applies to the zero temperature case only. 

Figure ^ shows F as a function of ^ea for p = 3 and several values of a as obtained 
from the numerical solution to Eq. ( |3.4.15| ). The T = transition takes place at 
the maximum value of F, F^. The corresponding value of ^ea is the discontinuity of 
the order parameter at the first-order transition. While F^ increases rapidly with 
a, the jump of the order parameter decreases as the strength of the coupling to the 
bath increases. The presence of the Ohmic bath thus tends to make the first order 
transition smoother. (This property tells us that it will be very difficult to see the 
first order transition by solving numerically the real-time dynamic equations.) 

This behavior results from the fact that (7i(e) diverges logarithmically as e ^ 0. 
In order to see this, we decompose the interval of integration in two parts, < x < 1 
and 1 < X < oo. The integral over the second interval is a finite constant at e = 0. 
In the integral over the first interval we may neglect x^ compared to x and write 

— (I-^ea)- r , = --lne + 0(l). (3.4.16) 

2 ^ ^^^^ .-0 Jo X + 2e + v^V^TIi 2 ^ ^ ^ ' 

We choose p = 3 for concreteness and solve Eq. (|3.4.16|) for F to obtain the equation 
of state in the high noise limit: 

1 

a' 



F = ^^3gEAe™(^-''-). (3.4.17) 
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Figure 5: The dynamic Edwards- Anderson parameter as a function of the couphng 
to an Ohmic bath, a, for F = 4 and (3 = 10. Sohd hue with point: exact numerical 
calculations. Dashed line: the low- frequency approximation of Eq. ( p. 4. 111) . 



This function has a maximum at 



where F reaches the value 



27ra 



(3.4.18) 



Tmax = ^ y exp ira . (3.4.19) 

We thus find the two features mentioned above, namely, a reduction of the disconti- 
nuity of the order parameter and a rapid increase of F^ for high values of a. Express- 
ing Eq. (|3.4.19|) in terms of the original variables of the problem (cf. Eq. ( p.l.7|) ) we 
find that, in the high noise limit, the T = dynamic freezing transition takes place 
at the critical coupling 

^^exp(-7ra) . (3.4.20) 

Thus, for a ^ 1, Jd is proportional to the exponentially small energy scale of 
Eq. (|2.21|) associated to incoherent tunneling in the isolated tls. It must be em- 
phasized that the existence of this scale is a feature of the spherical model used in 
this paper. Real tls {i.e., described by Ising spins) localize at a = 1. Therefore, 
is expected to vanish precisely at a = 1 for discrete spins. 

Deep in the ordered phase the system is expected to freeze with ^ea ~ 1- This 
regime occurs for sufficiently high values of a or sufficiently low values of F. Consider 
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Figure 6: F as a function of at T = for p = 3 and an Ohmic bath. The curves 
follow from the numerical solution to Eq. ( |3.4.15|) . The coupling to the bath a runs 
from to 1.2 in intervals of 0.2 from bottom to top. 



first the former case with Fa^ ^ 1. Then, e <^ 1 and we can still use Eq. ( |3.4.16| ) 
which, for ^ea ~ 1, reduces to 



^EA ~ 1 In 



Fa^ 



Ta^ > 1 



where Kp{l) = yp{p— l)/2. 

In the opposite case, Fa^ <^ 1, e is large. In this case 



(3.4.21) 



leading to 



gEA ~ 1 - 



F 



9i{e) 



1/2 



1 



3v/i V 8v/i 



+ 



3tc \Kp{l) 



F 



1/2 



a + 



(3.4.22) 



Fa^ < 1 . (3.4.23) 



In both regimes the effect of the noise leads to an increase in ^ea thus stabilizing 
the ordered phase. 



The sub Ohmic case 



Figure |^ shows F as a function of ^ea for p = 3 and several values of the coupling 
to a subOhmic bath with s = 1/2. The results were obtained by numerically solving 
Eq. ( p.4.11|) . The qualitative features of these curves are similar to those found in 
the Ohmic case. 

As discussed in Section ^, in the subOhmic case the isolated tls has a localization 
transition at a critical value a*^'^'^ of the coupling to the bath. We thus expect a 



20 



15.0 



10.0 




Figure 7: F as a function of ^ea at T = for p = 3 and a subOhmic bath, s = 1/2 
from Eq. ( 3.4.11| ). The couphng ai/2 runs from to 2 in intervals of 0.4 from bottom 
to top. 



transition to the ordered phase at J = for all a > a*^^'^ in the interacting system. 
Near the critical point at J = 0, e is small. For s < 1 the integral on the right-hand 
side of Eq. ( |3.4.12| ) is finite at e = and gs{0) = /s(0)/2 where fs is the function 
defined in Eq. ( |2.2UD . Detailed inspection of the behavior of gs{e) shows that, as 
e ^ 



e for < s < 1/2 

g^(0) - g,{e) oc { eln(l/6) for s = 1/2 

for 1/2 < s < 1 



(3.4.24) 



We consider for simplicity the case s < 1/2 and take p = 3 for concreteness. For 
a ~ a^"*, Eq. ( |3.4.12|) acquires the form 



JM 

IF' 



/gil oc 1 - 



a \ 2-s 



a' 



crit 



EA), 



(3.4.25) 



where a™'* is given in Eq. ( |2.22| ). There is a maximum at oc yJM/Ti^j . The 
dynamic transition thus takes place at 



J.oc- fl 



a 



a 



crit 



1/2 



(3.4.26) 



The jump of the order parameter at the transition is (1 — a/a™*). Therefore, 

for a = a'^'''* the dynamic transition is continuous. Summarizing, for a < a'^''" the 



transition between PM phase and SG phase occurs at a finite value of Jd while for 
a > a*^^'^ an infinitesimal J is enough to render the system glassy. We expect to 
obtain this same behavior for an interacting tls in an Ohmic bath. 
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At large couplings, a ^ a^^^*', Qea ~ 1 and we find 

(3.4.27) 



a 



1 

crit\ 2-i 



a 



Notice the absence of J-dependent corrections that appear at higher order {a ^). 
In the opposite limit (large J), e is large. Then, 

g^ie) = ^ (l - ^se^'^'-'^ + ■■■), (3.4.28) 

with a, a constant. We find 



4 

1 - 



1/2 



4 ^(1) V ^ / smf 



(3.4.29) 



As before, the presence of noise favors the ordered phase. The comparison of 
Eqs. ( p. 4. 211 ) and ( p. 4. 27] ) shows that, at strong coupling, an Ohmic bath is more 
effective than a subOhmic bath in freezing the spins. At weak coupling we have 
a linear dependence on a in both cases. For F <^ 1, however, the slope is higher 
in the latter case which results now in higher values of ^ea- Notice the presence of 

the extra factor (huph/J^ that amplifies this effect if the phonon energy is larger 
than the magnetic energy. 

The superOhmic case 

In the superOhmic case no localization transition exists at J = 0. For s > 1, 
5's(e) diverges as e~ in the limit e — >■ 0. This corresponds to small J or large as- 
A calculation similar to those performed above yields the critical coupling 

/sin— X''"^ 

J^r^hcupj—^j , a>l. (3.4.30) 

As in the Ohmic case, the critical coupling decreases with increasing a but only 
as a power law. The jump of the order parameter at the transition is however 
independent of a. 

Deep in the ordered phase, for small values of J, we find 

,_l_(^)^(!!lil)-'. ,3.,31) 

For small values of a Eq. ( ^.4.29| ) is still valid. Notice that for s > 1 the 
enhancement of the order parameter due to the coupling to the bath decreases when 
ojphl J increases. 
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3.5 The real-time correlation function 



In thermodynamic equilibrium the correlation function and the imaginary part of 
the susceptibility are related by 

1 doj 
C'(^) = T7 T.'^^i{^)s^m = qEA + h —x'Lci^) coth(/?^cu/2) cosM). (3.5.1) 

iV ~^ Jo TT 



If instead of the equilibrium response function we use in Eq. ( p.5.1|) the expres 



sion for x"{^) obtained through the AMS we obtain a correlation function that is 
closely related to the stationary part of that obtained through real-time dynamical 
calculations. This relationship was discussed extensively in CGS in the case of the 
isolated system with the following conclusions: 

1. The temperature below which the AMS solution exists coincides precisely 
with the dynamical critical temperature obtained from the dynamical calcu- 
lations. This is the temperature below which the real-time dynamics of the 
system becomes non-stationary and violations of FDT set in. 

2. The parameter m precisely coincides with X, the FDT violation factor. This 
is related to the effective temperature of the system in the aging regime, Tes = 
T/X H, see Section 4. 



The response function derived from the AMS is identical to the out-of-equilibrium 
response function when the long waiting-time is taken first and the weak- 
coupling limit in taken later on. More precisely, 

Cams(^) = ^y"('5i(t)si(0))AMs = lim lim CoYN{t + t^,t^) . (3.5.2) 



A proof of the analogous properties for the system coupled to the bath can be 
given following the same lines. The first two conclusions remain unchanged. The 
third one generalizes to 

CAMsit) = lim CoYNit + t^,t^) . (3.5.3) 

valid for all values of a. The aging regime, t t^, in which C^YNit+t^, t^) decreases 
below is not accessible in this approach. 

In this Section we shall analyze in detail several time regimes in C{t). We use 
throughout this section the original variables of the problem. 



3.5.1 No coupling to the bath 

We consider first the case in which there is no coupling to a bath. Then, the analytic 
continuation of Eq. ( |3.4.4| ) is 



X»„M = ^(jj) ^K- — . (3.5.4) 
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Substituting this expression in Eq. ( |3.5.1| ) and making the change of variables uj = 
4:KpJ/M X in the integrals we obtain the correlation function 

= ^EA + — f T7^) / dx X Vl- x^cos{xt/to) coth ( -z^x] , (3.5.5) 
TTKp \MJ J Jo I 



where to is a characteristic time given by 

to = . (3.5.6) 

At T = Eq. ( gXsD reduces to 



. 4:KpJ , 



C{t) = Oea H — ^ / dx X Vl — x'^ cos( Xt/tn) 

TlKp \MjJ Jo ^ ' ^ 

i^Y' iF,(l;l/2,5/2;-(t/to)V4) 



(3.5.7) 



^EA + 



where 1F2 is a generalized hypergeometric function. From the normalization condi- 
tion C{t) = 1 we find the quantum equation of state, 

K,(l-,EA) = 3^(^) . (3.5.8) 

found previously in CGS. The asymptotic behavior of the correlation function in the 
long-time limit is 

C^W - ^- + ;^(^)'^7^)'^'/(t/to), (3.5.9) 
t>to 37r/tp \MJ J \t J 

where / is an oscillatory function that can be expressed in terms of Fresnel integrals. 
From Eqs. ( p.5.8|) and ( p.5.6|) the frequency of the oscillations is 

^0 ~ jj-^ r . (3.5.10) 

M(l - ^ea) 

At the dynamic transition point ^ea depends only on p. Then, ujq oc Ti/M, the 
characteristic frequency of the non-interacting tls. Deep in the ordered phase ^ea ~ 



1 and Eq. (|3]|) yields 1 - ^ea ~ Ti/yMJ. Then, in this limit ujq ~ y J/M 



At temperatures higher than T^r = h/to ~ h^J J/M , but low so that the results 
from the approximation can still be used, we can approximate coth 2; ~ in the 
integral on the right-hand side of Eq. ( p. 5. 5 ) and write 



C{t) ~ gEA + / dx Vl — cos (— 

nJnp Jo \ to 



(3.5.11) 



Jun ytoJ t 
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where Ji (x) is the Bessel function. Notice that Eq. ( |3.5.11| ) also holds for all tem- 



peratures for times t^h/T. The normalization condition now yields 

Kp{l — ^ea) = — - I dx \J\ — = — , (3.5.12) 
vr J "'0 J 



which is the classical equation for ^ea F2|. In this classical regime the long-time 



asymptotic behavior of the correlation function is 

Notice that the power-law decay of the amplitude of the oscillations oc t^^/^ at high 
and low temperatures is the same. 

3.5.2 Finite coupling to a bath. 

In the presence of a coupling to an Ohmic bath there are two different regimes. 
At frequencies higher than = ha/M the inertial term in Eq. ( p.4.4| ) dominates 



over the term proportional to a. For times shorter than t* = l/ui* the system thus 
behaves as if it were isolated. At longer times, when inertia may be ignored, we 
have 



/ — / * \ 1/2 



naoo + ZKpJ/h \ zKpJ ) 
and the motion is overdamped. The correlation function then reads 

C{t) H — jpi \ du y/uj cos ujtcoth.{(3huj/2) , (3.5.15) 

71 ( 2KpJ ) 

where 70 = ah is the classical friction coefficient. Performing the integral we find 



2h^Jo/7T 



3 



-r2 , T = o 



C{t) - , X ^ (3-5.16) 

AKpJ j [ 4T/h t-2 , T > h/t . 

Notice the difference in sign between the results at zero and finite temperature. At 
zero temperature C (T) approaches ^ea from below, whereas at T 7^ it does so from 
above. In the Ohmic case the exponent controlling the decay of the T = correlation 
function is the same that controls the amplitude of the coherent oscillations found 
in the absence of noise. 

At finite temperature the decay is slower, C{t) — ^ea oc t"^/^ . In the classical 
model, the non-equilibrium correlation function C{t + tw, t^) approaches the plateau 

as C(t + t^, t^) — Qea oc t~'^^'^^ for t <^ t^. It was found that the temperature- 
dependent exponent z^(T) approaches 1/2 in the zero temperature limit in agreement 
with our result. The calculation of the temperature corrections to the exponent lies 
beyond the power of our low-temperature approximation. 
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At finite temperature, in the long-time limit, our results coincide with those 
obtained from the solution of the classical Langevin equation without inertia. Al- 
though the asymptotic form of the correlation function is independent of M {i.e., 
of the tunneling frequency A) it must be remembered that Eq. ( 3.5. 16| ) only holds 



for times longer than t* which does depend on A. A consequence of this fact is that 
the dynamics of the model in the limit A ^ is trivial. Indeed, it can be shown 
from Eq. (|A|) that 

XregH „ 7^ , jS-j ■ (3.5.17) 

Then, for any finite uj, 

hm XhecM^O. (3.5.18) 

However XregI"^) cannot be identically zero since the static susceptibility, Xreg(O), 
is finite and it is given by Xreg(O) = k"^ according to Eq. ( p.4.4|) . Xreg(O) can also 
be expressed as 

,,„(0) 24iM . (3,5,19) 

J-oo TT UO 



Equations ( |3.5.18| ) and ( |3.5.19| ) are compatible only if 



lim ^3^^^ = _ 8{uj). (3.5.20) 

M/7o->oo UJ Kp 

Therefore, the system has no intrinsic dynamics in this limit. In terms of the original 
spin model this is a simple consequence of the form of our starting Hamiltonian, 
Eq. ( p.4|) : if A = the spin variables commute with the Hamiltonian and are thus 
constants of the motion. In terms of the particle interpretation the limit M/ 70 oo 
corresponds to an infinitely massive particle that is not able to move or to the limit 
of zero friction where there is no dissipation. 

The expressions in Eq. ( |3.5.16| ) can be readily generalized to non-Ohmic baths. 
We find that in the long-time limit: 

r cos vr) , T = , 

C(t)-gEA oc I (3.5.21) 
[ cos (f tt) t-t , T > h/t . 

4 Real-time dynamics 

In this Section we study the real-time dynamics of the system coupled to the envi- 
ronment. We use the dynamic equations for the symmetrized correlation and linear 



response functions derived in |22| with the Schwinger-Keldysh formalism and we 
solve them numerically, as a function of time, for different couplings to the bath and 
different environments. We compare the results to the ones obtained in the previous 
Section with the imaginary time formalism. 
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4.1 The dynamic equations 

The dynamic equations for the model defined in Section 2 were derived in |^ . They 
are of the Schwinger-Dyson form and read 



{Mdf + z{t))R{t,t') 
{Mdl + z{t))C{tX) 



S{t -t')+ / dt" S(t, t")R{t", t') , (4.1) 
Jo 

roo rt' 

/ dt"^{t,t")C{t",t')+ dt" D{t,t")R{t',t") ,{A.2) 
Jo Jo 



with the equal-times conditions C(t,t) = 1 and R(t,t) = and 



lim dtR(t,t') = — 



Jiin dtR{t, t') 
lim dtC{t,t') 







lim dtC{t, t') = Q. 



(4.3) 
(4.4) 



The equation for the Lagrange multiplier, z(t), reads 



z{t) 



dt" [S(t, t")C{t, t") + D{t, t")R{t, t")] 

+M f dt" f dt"' {dtR{t, t")) D{t", t"') {dtR{t, t"')) 
Jo Jo 

+ \dtR{t, s)dlC{s, t) - dlR{t, s)dt'C{s, t' 



(4.5) 



s ^ 
t^t' 



The total self-energy and vertex include the interaction with the bath and are given 
by 



S(t,t') = -Ar]{t-t')-^\m 
D{t,t') = 2hi^{t - t') + ^Re 



-lp-1 



C(t,t')-yi?(t,t') 
C{t,t!)-%{R{t,t!)^R{t!,t)) 



with 



v{t - t') = dujl{uj) coth cos(cu(t - t')) , 

7](t-t') = -Oit-t') / duj liu) siniujit-t')) . 

Jo 

The spectral density of the bath, I{uj), has been defined in Eq. (|2.11| ) p8 



(4.6) 
(4.7) 

(4.8) 
(4.9) 



In the following we shall compare the effect of environments with different val- 
ues of s and using different coupling strengths. The high-frequency cut-off, Uc, is 
introduced to avoid the divergence of i^{t). In the Sub-Ohmic case, when we solve 
the equations numerically, we also need a low frequency cut-off, that we impose in 
a hard manner by including a factor 9{uj — h) in the definition of /(t^). 
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The kernels v and r/ can be computed for all values of s. In the numerical solution 
to the set of coupled integro-differential equations ([4.1|)-( ^72D it is more useful to use 
the integral of the kernel 77, //(r) = J^dT'r]{T'), that reads 

and, when s takes the values 1/2, 1, 3/2, it becomes 

s = 1/2 subOhmic 

s = 1 Ohmic 
s = 3/2 superOhmic 

On the right-hand-side we have written the limiting form for uJcT ^1. It is clear 
that, as for the imaginary-time kernels, the dependence on Uc is very different in 
each of these cases. 

We shall rescale the real time and the other parameters and functions in the dy- 
namic equations to match the definitions that we used in the Matsubara calculation. 
Under the rescaling of time, t — ^ J/ht, the correlation function remains unchanged 
and the response transforms as R ^ hR. The rescaled dynamic equations are 
identical to Eqs. (|4.1| ) and ([4.2|) with M replaced by F^^. 




4.2 Numerical study of the real-time dynamics 

As shown in Section |3.2| , both static and dynamic transition lines depend strongly 
on the strength of the coupling between system and bath. We can also see this effect 
by following the real-time dynamics of the system coupled to the environment. We 
have solved Eqs. ( [4.1|) -( ^l2D numerically with a predictor-corrector algorithm that 
allows us to reach long times with a high accuracy. For each set of parameters 
we have checked the data collapse for different values of the iteration step h in the 
discretized equations. In general, there is a good collapse for h < 0.02 and, typically, 
we have used h = 0.01 and h = 0.02. 



4.2.1 Effect of the interactions: localization against glassy behavior 

In the Introduction and Section 1 we recalled several results for localization in dilute 
two-level systems coupled to a bath. In this paper we focus on a soft spin version 
of the interacting problem. Our first aim is to determine the effect of the coupling 
J on the localization properties of this system from a real-time dynamic point of 
view. In Fig. |^ we show the decay of the symmetrized correlation C(t + t^, t^) using 
a subOhmic bath with s = 0.5, Uc = 10 and Uph = 5. The three upper curves were 
obtained for a = 4 and changing the value of the SG coupling strength J. When 
J = the system localizes for a > a*^^^^: for any ty^ and long enough t the correlation 
reaches a plateau and it does not decay below this value. When a small coupling 
is switched on the decay changes. The correlation approaches a plateau for small 
values of t — tu, but it subsequently quits the plateau and decays towards zero. The 
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Figure 8: The decay of the symmetrized correlation in three cases: localization for 
a = 4 and J = 0, glassy decay for two nonvanishing values of J, J = 0.5 and J = 1 
and a simple decay towards zero for the case of a small coupling to the bath, a = 0.2 
and J = 0. We have chosen a subOhmic bath with s = 0.5, Uph = 5 and Uc = 10. 
The quantum parameter F equals one. 

system has glassy non equilibrium dynamics that we shall quantify below. Finally, 
when the coupling to the bath is very small and J = the system does not localize 
and the correlation decays to zero with wide oscillations. 

4.2.2 Dynamics in the paramagnetic phase 

For a chosen coupling to a bath, at sufficiently high values of F and/or T the system 
equilibrates with the environment and it quickly reaches a stationary regime where 
the quantum fluctuation - dissipation theorem (fdt) is satisfied. This property has 
been proven for the p spin model in [^|, for the large N fully connected Heisenberg 



SU(N) model in and for a soft version of the quantum model in |^9|. In all cases 



the systems were coupled to an Ohmic environment and the limit of weak coupling, 
lim„^o linii„->oo, was considered. The correlation and response have a rapid decay 
towards zero with oscillations that depend on the value of the quantum parameter 
F and, as we show here, on the coupling to and the type of environment used. 

In this Section we analyze the effect of the quantum fluctuations and the bath 
on the conclusions mentioned above. We flrst consider a flxed Ohmic environment, 
i.e. we take s = 1 and we flx Uc = 10. We display in Fig. R| the decay of the 
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Figure 9: The stationary autocorrelation (left) and response (right) functions at 
T = 2 for several values of the quantum parameter F, given in the key. The bath is 
Ohmic and lJc = 5, a = 0.8. 



correlation and response functions for different values of F in the PM phase. It is 
clear from the figure that the period of the oscillations decreases with F. In order 
to quantify this dependence one can compute x"{^)/^ and follow the evolution of 



the peaks. We show two examples in Fig. 10. The data on the left panel correspond 



to those on Fig. |^. On the right panel we represent x"(a;)/a; for T = 0, F = 5 
and several values of a. For small values of a the system is deep in the PM phase 
and there is a well defined peak in x"(u;)/u; at a finite value uo that increases with 
decreasing a. At high enough values of a a tail at low frequencies starts developing, 
indicating that the dynamics is slower and that the system approaches the transition 
towards the glassy phase. Eventually, as discussed in Section P?T| , for high enough 
a, the parameters fall below the transition and the system becomes glassy with slow 
dynamics. 

4.2.3 Dynamics in the glassy phase 

In Figs. we compare the behavior of the correlation and response functions in the 
glassy phase when the system is coupled to an Ohmic environment through different 
coupling constants. We choose T = 0.1, F = 4 and we compare the effect of a = 0.2 
and a = 1. The high-frequency cut-off is Uc = 5. From the discussion in Section 3 
we expect that the system is in the PM phase in the first case and in the SG phase 



in the second. This is seen in Fig. |TT]. For a = 0.2 the correlations rapidly reach a 
stationary regime and they oscillate around zero. For a = 1 the behavior is different. 
There is a first rapid decay towards a plateau, that has a low value, and then, a 
slow and monotonic decay towards zero. Aging effects are apparent from the figure. 
The response function also shows a qualitatively different behavior according to the 
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Figure 10: The frequency dependence of x"{l^)/^- Left panel: T = 0, a = 0.6, 
= 5 and several values of F. Right panel: F = 5, T = 0, = 10 and several 
values of a. The environment is Ohmic in both cases. 



value of a. In one case it quickly acquires a stationary oscillatory behavior around 
zero, in the other it has a long tail as expected in a glassy system. We then conclude 
that the system has undergone a dynamic phase transition between the PM and SG 
phases at an intermediate value of a. 

An approximate expression for the dependence of the Edwards- Anderson param- 
eter on F and a, at T ~ 0, has been obtained in Section 3^. We can also check 
this law by estimating the value of from the numerical solution of the real-time 
equations. If we plot the correlation function for several values of a in a log-log 
scale the plateau at ^ea can be easily identified. It is a slowly growing function of a 
that is rather well described by Eq. ( |3.4.21| ). 



We also investigated the effect of different environments (different s) of the same 
strength (same a) using the same value of the high-frequency cut-off that we took 
equal to ujph- From the discussion in Section 3.3 for some values of ujph we expect 
the relaxation to be slowest for the subOhmic bath, intermediate in the Ohmic case 



and faster for a superOhmic environment. This is illustrated in Fig. 12. The decay 



is slower when s = 0.5 than in the other cases. In the extreme case of s = 4 the 
system has gone across the transition towards the PM phase. However, this behavior 
is not generic. 

The relation between the correlation and response plays a key role in the descrip- 
tion of the dynamic behavior of glassy systems. When the system is in equilibrium, 
this relation is model independent and it is given by the fluctuation-dissipation 
theorem (fdt). When the system is glassy and it evolves out of equilibrium, the 
conditions to prove the theorem are not satisfied but simple generalizations have 



been exhibited in a number of models |42, 22 
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Figure 11: Left: The auto-correlation C{t + t^j^tw) as a function of t for a = 0.2 
(PM phase) and a = 1 (SG phase). The temperature is T = 0.1 in both cases. The 
different curves correspond to different waiting times, tyj = 5, 10, 20. For a = 0.2 
the curves collapse on an asymptotic one, while for a = 1 they show aging effects. 
Right: The response R{t + t^,t^) as a function of t for the same parameters. The 
effect is similar. 




Figure 12: The symmetrized autocorrelation and the linear response at T = 0.1, for 
r = 1. We compare the effect of a sub-Ohmic (s = 0.5), an Ohmic (s = 1) and a 
two super-Ohmic (s = 1.5, s = 4) baths. See the key for the details. The coupling 
to the bath is kept fixed to a = 3 and the high-frequency cut-off equals the phonon 
frequency, Uc = ojph = 5. 
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The quantum FDT for a system in equilibrium, in the rescaled variables, reads: 



R{t) = e{t) I p duj e-'"* tanh [^-^ C{u) (4.11) 

where _ 

C{Co) = 2Re / dr e*'^* C(r) . (4.12) 
Jo 

The quantum FDT is an integral relation between the stationary linear response and 
the symmetrized correlation function. 

The asymptotic dynamics in the glassy phase take place in two time scales that 
are separated by the plateau in the correlation function. As shown in for the 



weak coupling limit, the stationary part of the decay, when the correlation decays 
from 1 to is such that the quantum FDT holds. We have checked that this result 
also holds when the system is strongly coupled to a nonOhmic bath. In the weak 
coupling limit, when the correlation decays beyond EA, the relation between linear 
response and correlations takes the form of the classical FDT and it reads 

R{t) = 9{t) (3,,J ^C{t) . (4.13) 



where Peff is the inverse of an effective temperature |Q TgFp and J appears since 
we have rescaled time. A concrete way of testing the validity of this equation is to 
plot the integrated response function 

ft+tw 

X{t + t^,t^)= dt'R{t + t^,t') (4.14) 



against the symmetrized correlation C{t + tw,tw) for a long using t as 

a parameter. For short time-differences, when t— ^ and Cit+t^, t^) > Qea, this 
construction does not have any particular meaning and the curve is nonmonotonic 
with strong oscillations. Instead, when t — tw ^ ty^ or longer and C{t + tw, tw) < Qea, 
the curve becomes a straight line of slope — 1/Teff- 

In Figs. |1^ we display the x against C plots for different values of the parameters, 
explained in the caption and keys. The panel on the left shows the x C curve for 
different values of the coupling a. The slopes of the curves, and hence Teff, change 
smoothly and there is a clear non trivial dependence on this parameter. The panel 
on the right displays the x{C) plot for a fixed value of the coupling a and several 
values of s, s = 1,1.5,4. The first two cases are in the glassy phase while the latter 
falls in the PM phase and the parametric plot does not show a straight line piece. It 
is difficult to decide from these figures if the slopes depend on s or not. 

In order to sharpen our conclusions about the dependence of Teff on the char- 
acteristics of the environment we take profit of the empiric relation between Teff 
and the breaking point parameter m in the replica analysis of the same model, 
Teff = T/m. In Section 3 we developed a low temperature, low frequency approxi- 
mation to solve the saddle point equations stemming from the replicated Matsubara 
analysis of this problem. In these limits we derived a set of equations that link 
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Figure 13: Left: The parametric x against C curves for different values of the 
couphng to the environment, at fixed T = 0.1, tUc = 5 and using an Ohmic bath. 
Right: the dependence of the x against C plot on the kind of bath used. We include 
a straight line as a guide to the eye. In all curves T = 0.1, a = 2 and c^c = 5 for all 
curves. 

T/m to a and s that can be solved numerically. We found that for fixed s the 
effective temperature Teff is a growing function of a. This result is reminiscent to 
the dependence of Teff on the external temperature T in a classical problem: the 
lower T, the higher Teff meaning that higher values of the effective temperature 
are reached when the system is deeper in the glassy phase. Each curve approaches 
one when a oo and the corrections can be read from the asymptotic analysis 
presented in Section 3. The dependence of Teff on s is weak but non-monotonic. 
(We have already encountered a non-monotonicity related to the fact that the factor 
ujph^ changes the coupling between system and bath differently for different values 
of s.) 

5 Conclusions 

In this article we discussed the effect of a quantum environment on the nonequilib- 
rium dynamic properties of an interacting quantum glassy system. We have shown 
that, as in the case of a simple tls, the influence of the quantum bath is very 
important. 

Two limits of the quantum model are easy to derive or were already known. On 
the one hand, in the absence of interactions, the calculations shown in Section 2 
and the numerical results of Section 4.1.1 prove that when the model is coupled to a 
subOhmic bath, it undergoes a localization transition at a critical value of the cou- 
pling a. The localized phase is characterized by a symmetrized two-time correlation 
function that, as a function of time-difference, approaches a non-vanishing asymp- 
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totic value and never decays to zero. On the other hand, it was known that when 
interactions are switched on and the hmit of weak couphng, a — > 0, is taken, the 
model has glassy dynamics with a symmetrized correlation function that depends 
on the waiting-time and decays in two-steps with a first approach to a plateau and 
a second decay towards zero p2| . 

The aim of this article was to analyze the combined effects of the interaction 
(J 7^ 0) and a strong coupling (a 7^ 0) to quantum environments of different types 
(different s). We summarize our findings as follows: 

First, we determined if the model has a localized phase in the presence of inter- 
actions. How to define such a phase for an interacting system is a difficult question 



(see, e.g., [^). Here, we adopted as evidence for a localized phase the fact that 
for a long enough waiting-time ty^ the correlation function does not decay to zero at 
any time-difference t — t^. With this criterium we saw that, as expected, there is no 
localized phase when interactions are switched on. 

This result can be interpreted a posteriori by resorting to the concept of effective 
temperatures generated by the nonequilibrum dynamics of glassy systems. Indeed, 
it has been shown for classical systems that the modification of the fluctuation- 
dissipation theorem observed in systems evolving slowly out of equilibrium is related 
to the self-generation of effective temperatures (typically higher than the one of the 
environment) The proof presented in has not been extended to quantum 
systems yet. However, as argued in |^ for the quantum model studied in this paper 
when weakly coupled to an environment, the slow part of the relaxation looks clas- 
sical, with a quantum fluctuation - dissipation relation that became classical with 
an effective temperature that is higher than the temperature of the environment. 
In particular, when the model is coupled a quantum bath at zero temperature it 
acquires a non-vanishing effective temperature. This effect has been observed in 
other quantum glassy systems too ||2^, When the system is strongly coupled 
to the environment the relaxation slows down with respect to the weakly coupled 
case. However, the two step relaxation remains with a slow regime controlled by 
a nonvanishing effective temperature. Thus, we conclude that the generation of an 
effective temperature by the interactions is consistent with the fact that the system 
does not localize. It is well known that even in simple tls the localization effects 
disappear at finite temperature. 

Next , we analyzed the effect of a strong coupling to an environment on the glassy 
properties of the model. We showed that stronger couplings to the bath favor the 
glassy phase for any type of bath. By this we mean that for larger value of a the 
area of the spin-glass phase on the (T, F) plane increases. We also characterized 
the dependence on a of several properties of the system as the Edwards- Anderson 
parameter, the effective temperature, etc. 

Finally, we studied the effect of different types of baths. Concerning this issue 
the conclusions are cumbersome given the fact that a new parameter, the phonon 
frequency Uph, appears in the spectral density when s 7^ 1. If Uph is not equal to one, 
the effect of different baths are complicated. For instance, the dependence of ^ea on 
s can be nonmonotonic as well as the location of the critical line on the (T, F) plane. 
We exhibited some examples but we cannot draw general conclusions concerning 
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this issue. 

The motivation for this study were manifold. The effect of quantum environ- 
ments on interacting macroscopic quantum systems is a problem that is now being 
revisited in the context of quantum computing ||0]. Decoherence, or how quantum 
interference effects are lost due to the interaction with the environment, has to be 
as much reduced as possible to make a quantum computer performant. Again in the 
context of quantum computing, an isolated Edwards- Anderson quantum model in a 
random transverse field has been proposed to mimic an isolated quantum computer 
with (short-range) interactions between the spins (that represent qubits) and with 
static "imperfections" in the individual two-level system energies In this work 
we analyzed a soft limit of a disordered quantum model with long-range p > 2 inter- 
actions in a transverse field. It would be very interesting to see which, if any, of our 
conclusions are modified if the soft spin limit is lifted and, even more importantly, 
if a finite dimensional model is considered. This project, however, is a very difficult 
one. 

On a more physical side, glassy phases at very low temperatures where quantum 
fluctuations are important have been identified in a number of physical systems. In 
the proper analysis of these systems the role played by the quantum environment has 
to be taken into account. Our results are a first step towards the characterization 
of the effects of the environment. Again, it would be interesting to go beyond the 
mean-field limit and derive similar results for a finite dimensional model. 
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